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The correlations in the spectra of quantum systems are intimately related 
to correlations which are of genuine classical origin, and which appear in the 
spectra of actions of the classical periodic orbits of the corresponding classical 
systems. We review this duality and the semiclassical theory which brings it 
about. The conjecture that the quantum spectral statistics are described in 
terms of Random Matrix Theory, leads to the proposition that the classical 
two-point correlation function is given also in terms of a universal function. 
We study in detail the spectrum of actions of the Baker map, and use it to 
^ \ illustrate the steps needed to reveal the classical correlations, their origin and 

£^ , their relation to symbolic dynamics. 
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^ ! 1 Introduction 

O . 

The interest in the correlation between classical periodic orbits, and in particular, in 
the spectrum of their actions, emerges from the attempts to provide a semiclassical 
proof of the universality of spectral correlations of quantum chaotic systems ( "the 
BGS conjecture" Action correlations were first discussed by Argaman et, al. 
. | were the universality of action correlations, and their relation to Random Matrix 

^ . Theory (RMT) were studied for a few chaotic systems Various aspects of the 

subject were investigated later [[| [|, |[ [?], [|, ||, Q. This culminated recently in the 
work of Sieber and Richter (ll[] , who identified pairs of correlated trajectories, whose 
contribution to the spectral form factor for systems with time reversal symmetry 
is identical to the next to leading order term in the formfactor predicted by RMT 
(see also @,|l|). 

The purpose of this paper is twofold: first, to review the semiclassical context 
where action correlations and their expected universal features arise in a natural 
way H This will be done in the main body of the present section, where 

also the connection with RMT will be made explicit. Second, to test the general 
semiclassical arguments on action correlations for a paradigm chaotic dynamical 
system - the Baker map. This system was investigated previously by a number of 
groups, (2[ H| |6[ |8), who demonstrated numerically the existence of the expected 
correlations. Here, we develop another approach for the analysis of the action 
spectrum, where we try to systematically asses the way the periodic orbits and 
their actions can be partitioned to families which are dynamically related. One 
aspect of this approach is studied by casting the problem as an Ising model (in one 
dimension and with a long range, yet exponentially decaying interaction). Moreover, 
we show several features of the action correlations which escaped the attention 
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of previous works, and in particular, we analyze the correlations in terms of the 
symbolic dynamics. We hope that this insight will pave the way to a more complete 
understanding of the universal features of action correlations in general. 



1.1 Action correlations and the semiclassical theory of spec- 
tral statistics 

Consider a finite, two dimensional domain A and an area preserving map 

M : 7' =M(<y) ; j',je A. (1) 

Any area preserving map can be expressed implicitly through a generating function 
|T^| . Let 7 = (x,y) , 7' = (x',y r ). The generating function Fi(x',x) defines the 
map through the relations 

dx dx' 

where x' is to be expressed in terms of (x, y) by solving the first equation, and y' 
is given explicitly using the second equation. Twist maps are maps for which the 
first equation above has a unique solution for any (x, y) = 7 G A. From now on we 
shall assume M. to be a hyperbolic twist map. 

Periodic points of period n are solutions of the equation 7 = A4 n (-f). The 
corresponding rt-periodic orbit is obtained by iterating the map: 

The action assigned to a periodic orbit is 

n 

f a = J2Fr(4 a) > x ^) ; 4 8) = *£° (3) 
3=1 

The number of n-periodic orbits, N n increases exponentially with n. The object of 
our investigations is the set of actions {f a }a=ii in the limit n — > 00. 

Under quite general conditions one can show that the mean action (/) and the 
variance var/(n) are both proportional to n. Moreover, the actions of n-periodic 
orbits are bounded within an interval which grows algebraically with n. However, 
since their number grows exponentially with n, one expects exponentially small 
spacings between successive actions. We shall explain now why the pair correlations 
of action are expected to be universal, and determined by RMT. To this end we 
should address the quantum analogue of the map, in the semiclassical limit. 

The quantum analogue of the classical map, is a unitary evolution operator Un 
which acts on a N dimensional Hilbert space, with 



N 



2nh 



(4) 



where [•] stands for the integer value and | \A\ \ is the area of A. 

The spectrum of Un consists of N unimodular complex numbers e l6 '( N \ 6i(N) € 
[0, 2?r). Ample numerical evidence supports the conjecture that the spectral statis- 
tics of Un, is well reproduced by the predictions of random matrix theory (RMT). 
The quantum spectral density is denoted by 



d qm (0;N) =^2s{e-9i(N)). (5) 



1=1 
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In the semiclassical approximation, 



(x'\U N \x) 



2nhi 



d 2 F 1 (x',x) \ 2 iFl[x ' tX)/h 



dx'dx 



(6) 



Here, x is a discrete index which labels the eigenstates of the position operator and 
Fi(x\ x) is the classical generating function. A detailed discussion of the semiclas- 
sical evolution operator and its properties can be found in j^] . 
The semiclassical approximation (@) yields 



&6 e me d qm (6; N) = trll% » V A% 
Jo aev n 



) e ifa/H+i^r a iJ, a 



where V n is the set of n-periodic orbits. are the stability amplitudes 

n 



A (n) = 



r a \det(I~T^ 



(7) 



(8) 



Here, r a stands for the repetition number if the periodic orbit is a repetition of a 
primitive n a -periodic orbit (n a — ^-). T a is the monodromy matrix and [i a is the 
Maslov index for the primitive orbit. 

At this point we introduce the classical density of actions of n-periodic orbits, 



dd(f;n) = 4 n) e ifr "^«(/-/«(n)), 

a6P„ 



and we obtain 



trf/£ 



df e l f/ h dd(J;n) 



(9) 



(10) 



Using (||), and measuring the actions in units of the phase space area s = f/\\A\\, 
we get the relation between the quantum and the classical densities: 
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AO e lu6 d qm {6-N) 



ds e lks d c i(s; n) 



(11) 



k=27rN 



This equation expresses the quantum - classical duality. It relates two densities 
which are very different: The quantum density gives a unit weight to all the eigen- 
phase on the unit circle, while in the classical density, the actions are weighted by 
the stability amplitudes, and are assigned with a "charge" ±1, ±z depending on 
the Maslov index. The duality relation is expressed via a Fourier transform which 
involves both the variables and the parameters which specify the quantum and the 
classical distributions. On the quantum side, N is a parameter which specifies the 
value of h, while n is the value of the variable conjugate to 8. On the classical 
side, n is a parameter which specifies the period of the ensemble of orbits under 
consideration, while N determines the value of k - the variable conjugate to s. In 
the sequel we shall always reserve the first position to the spectral variable or its 
conjugate, while the second position is reserved to the parameter which specifics 
the system. 

We shall focus our attention on the spectral formfactors, which are the Fourier 
transforms of the pair correlation functions. The quantum and the classical form- 
factors are given explicitly as 



K qm {u-N) 



1 

N 



N 

E 



i.i' 



1 

N 



N 

E 



^2TTiT{ei-9 v )d N 



(12) 
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where, r = jfe and = ^- is the mean density and v must be an integer since d, 
is a periodic function. The classical formfactor is, 

K c i{k;n)= J2 4" ) 4" )eif(r '"'^"'"' )ett(s '" v) - ( 13 ) 

a, a' — 1 



It follows from (11) that the quantum and the classical 2-point formfactors are 
related by 

K qm {v = n;N)^ ^K cl (k = 2nN;n) . (14) 

Hence, the quantum spectral correlations are reflected in the classical spectrum, 
and vice versa. This equation expresses the important semiclassical result that the 
quantum formfactor is obtained from the classical one by interrogating the spectrum 
of action differences on the scale of iV -1 ~ %. 

Equation ( |l4|) has to be understood in the sense of distributions, since the form- 
factors, as defined above, are the Fourier transforms of a sum of S functions. Indeed, 
the formfactors computed for a given system, fluctuate, and do not converge to a 
limit when N — ► oo. The appropriate way to overcome this difficulty in the present 
context is to apply a smoothing procedure, which enables the extraction of well 
defined limit distributions. 

Starting with the quantum formfactor, we adopt spectral averaging which is 
based on the assumption of "spectral ergodicity". We order the spectrum so that 
8i < Oi+i, and partition it to N g subsets a g , each consisting of N — subsequent 
phases. Neglecting correlations between phases in different subsets, we rewrite 
as 

K qm {r;N) * ± £ I £ j^-O,)^ 

9 g=l iV l,l'<£cr g 

= ^E^( r ^)- ( 15 ) 

9 3=1 

K q m(T; N) is the formfactor for the spectrum obtained by multiplying all 9i € a g by 
^ so that they cover uniformly the entire circle. Taking now the limit N — > oo at a 



constant r, with N g » N rj v N, it is expected that K qm {r\ N), which is expressed 
now as an average over the ensemble of the partial formfactors, converges to a limit 
distribution K qm (r) which reproduces the prediction of RMT for the appropriate 
ensemble. This procedure is justified by the fact that in the quantum spectrum, the 
correlation range is of the order of a mean spectral spacing. Hence, the elimination 
of the correlations between different sets in ( |i~5j ) introduces a small error, which 
vanishes in the large N limit. 

A different smoothing procedure is required for the discussion of the classical 
spectrum. As will be shown in the sequel, the classical correlation length is of order 
of the inverse of n which exceeds by far the mean spacing which is exponentially 
small in n. A spectral smoothing should therefore rely on a different partitioning of 
the period orbits, such that the relevant correlations are preserved within a subset, 
while members of different subsets are statistically uncorrelated. One of the main 
problems in dealing with the classical correlations is the proper definition of such 
subsets 0. The smoothing of the classical spectra is by far more important since 
the number of actions increases exponentially with n. 

If we denote the results of the quantum and the classical smoothing procedures 
by (•), we obtain from ( |i~I| ) the relation which is basic to the present approach, 

1 /•,, n s> n 
— (K cl (2n-\n)) ; t= — 
N t N 



(K qm (T;N))^-(K cl (2n-;n)) ; r = -. (16) 
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Assuming that (K qm (n;N)) follows the RMT predictions, namely, (K qm (n; N)) — 
Krmt( t — n /N), we find that the classical spectra of chaotic systems must display 
universal pair correlations which can be derived from RMT using the relation ( |l6| ) . 
In particular, since for large TV, (K qm (T; TV)) is a function of r only, we derive two 
immediate predictions for (K c i(2tt— ; n)) : 



(i) (K c i(2n-] n)) is proportional to n, since the rhs of (16) must depend on TV only 
through the ratio jr. It is convenient to define 

K n (k) = ~(K cl (k;n)) , (17) 
n 

where the proportionality of (K c i(2tt^;ti)) to n is made explicit. 

(ii) /C n (fc) depends on k only through the scaled variable ^, 

tC / 77" 11 

JC n (k) = — K RMT ( — ) . (18) 
Zirn k 

Since k is the parameter conjugate to the action differences, it follows that the 
correlation length of the action spectrum is proportional to —. Moreover, ( |l8| ) 
shows that the classical correlations of actions corresponding to different periods 
are identical up to a scaling, and universal, since they are expressed in terms of 
a system independent function. These predictions pertain exclusively to classical 
properties and, once they are derived within classical mechanics, the road would be 
cleared for a semiclassical derivation of the BGS conjecture. 

The validity of (i) in the limit r — > was first shown by M. Berry |20). In the 
present context, it follows from the fact that the actions are discrete variables, and 
therefore, for sufficiently large TV the only correlations are the diagonal ones. Hence, 
for k — ► oo, 

(K cl (k;n)) - 5> a 4")| 2 « n(g) (19) 

a 

where g a is the number of different orbits which share the same action (mostly due 
to symmetry) and the summation is now restricted to trajectories with different 
actions. The equality on the right hand side is derived by using the ergodic sum- 
rule |Q and denoting by (g) the average value of the g a . Since (g) takes the value 1 
for systems without time reversal symmetry, and 2 for systems which are symmetric 
under time reversal, the leading terms in the RMT results for the CUE and COE 
are indeed reproduced. 

To derive the quantum formfactor in a consistent semiclassical way, one has to 
consider the classical formfactor for a fixed n, and change k (or TV) so that an interval 
of r values is scanned. However, to have meaningful pair correlations, TV should be 
larger than 2, and therefore the range of allowed r values, is at most < r < n/2. 
A further restriction on the range of r is due to the limited semiclassical accuracy, 
which is of the order of a mean spacing, |pl[ . In other words, the semiclassical trace 
formula may have its poles away from the unit circle, thus replacing the sharp S 
functions in (||), by spikes with a finite width, of the order of the mean spacing. 
Therefore, the regime r > 1 is not expected to be accessible to the semiclassical 
approximation, and thus, the range of applicability reduces further to the domain 
n < TV. It is important to note at this point that the required order of operations 
(keeping n fixed and taking TV — > oo) is consistent with the correspondence principle 
and the rules of quantum mechanics where the limit ?i — ^ is taken at a fixed value 
of the classical parameters. 

In the present work we shall study the Baker map, which will be introduced in 
the next section. It is a simple dynamical system, for which various properties can 
be derived analytically, yet it carries the full complexity of the action spectrum, and 
therefore it is appropriate as a paradigm. In section (||) the spectrum of the actions 



■5 



will be discussed from various points of view, and on the different relevant scales. 
The expression of the actions in terms of the symbolic codes of the periodic orbits 
will play a central role. Mapping the computation of the spectral density onto an 
Ising model enables us to introduce a few approximations which enable the deriva- 
tion of some properties of the action density ( |3.3| ) . The classical pair correlations 
are discussed in section (^). We study first the m'th neighbor distributions. We 
show numerically that as long as m is smaller than the action correlation length 1 /n 
measured in units of the mean action spacing, the m'th neighbor distributions are 
essentially Poissonian. Turning to the formfactor, we discuss alternative smoothing 
methods and relate them to the symbolic codes. We finally show that the classical 
correlations reproduce the expected RMT behavior as was conjectured in the gen- 
eral discussion presented above. We conclude the paper by a summary, where the 
connection between the present and previous works is discussed. 



2 The Baker map 
2.1 The mapping 

The Baker map is one of the most simple examples of chaotic maps [0, fl8| . 
It is an area preserving map of the unit square onto itself defined by 

x' = 2x-[2x] ; y'=\y+^[2x] (20) 

The stretching in the x direction and the squeezing in the y direction are responsible 
for the hyperbolic character of the map, while the "cutting and putting on top" gives 
the mixing property. Figure 0. shows one iteration of the map. 



Figure 1: The baker map 



The action of the map is easily translated to a Bernoulli shift: Every phase 
space point is presented in a binary basis 



y 



= ^^cii2 1 = 0-010203. 

i=l 
00 

- ]>>2- 1 =0- 6x6263.. 



i=l 



with Oi, hi S {0, 1}. The dynamics is given by shifting the binary point to the right 
when the two fractions are put back to back 



.636261 • 010203- 
x = • 010203. 
y = • 616263. 



...6 3 6 2 6iOi • 0203. 
x = • 0203... 
y' = • oi6i6 2 6 3 .. 



G 



2.2 Periodic orbits and codes 



A 71-periodic point is represented as an infinite repetition of a finite binary string 
with n entries v — {1/1,1/2, ....u n ),Ui € {0, f}. The cyclic permutations of (yi....v n ) 
represent the periodic points which constitute the periodic orbit. The number of 
n-periodic orbits is ~ 2™/n. The phase space coordinates of a n-periodic point can 
be written in term of the code 

^ n n 

■'■ t—ttY. 1 '--. 2 ' ; «• t—ttY,"-- - " ' : (21) 

3=1 3=1 

As n — ► oo periodic points fill phase space densely and uniformly. 

The symbolic dynamics introduced above is based on the partition of the unit 
square into two equal rectangles along the line x = 5. The sequence of binary 
symbols indicates the order by which the orbit visits the rectangles. Alternative 
codes can be generated by partitioning the unit square along the lines Xj = 
with r > 1 integer, and 1 < j < 2 r — 1. A symbolic code consisting of the 2 r symbols 
0, 1, • • • , 2 r — 1 indicates the order by which the periodic orbit visits the rectangles. 
The translation of a binary sequence to the 2 r code is done by considering successive 
sequences of r binary symbols as integers in [0, 1, • • • , 2 r — 1]. (Example, the binary 
code {001110101} is translated to the 2 r = 3 code {137652524}). Increasing r the 
code becomes more informative, because it locates the rectangles with an accuracy 
2 _r along the x axis. However, this is achieved at a cost: not all sequences of 
symbols are allowed. They are restricted by a Markovian grammar with a 2 r x 2 r 
connectivity matrix with only two non vanishing entries per row. 



C n, = \ 1 fOT (i'-2j)mod2-e{0,l} 
j j I otherwise 



The refined codes will be used in the sequel. 



2.3 Symmetries 

The mapping possess two discrete symmetries |l6| ]. The first is a space reflection 
symmetry 

R : x — > 1 — x and y — > 1 — y 

geometrically it is a double reflection about both the y = 1 — x diagonal and the 
y = x diagonal, it manifests itself on the code by R(vi) = 1 — Pi. The second is 
time reversal 

T : x ► y ; y > x and t — ► —t 

where reversing the time means reversing the mapping. Geometrically it is a re- 
flection about the y — x diagonal and its action on the code of a periodic orbit is 
T{vi) — Vn—i+x. Figure |[ shows the action of the symmetry transformations: R, 
f and RT = TR on the code v= (0 1 1 1). 



2.4 Generating function, action 

The action s associated with a phase space point (x, y) is defined as the generating 
function F\(x, x') of the mapping, with x' = x'(x,y). Because of the fact that in 
the Baker map, the x dynamics is independent of y, one has to derive the action by 
first extracting the generating function F2{x,y') from the conditions: 

dF 2 (x lV ') = t dF 2 (x 7 y') = ^ 

dx dy' 
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and get F2(x,y') = 2xy' — x[2x] — y'[2x] + const. To obtain F\{x,x') a Legendre 
transform is needed (see ]24|) 

Fi(x,x') = F 2 {x,y') -x'y' = -x[2x] . 

One should emphasize that F\(x, x') is not a generating function of the mapping. 
However choosing F\(x, x') as the action is consistent with the action obtained from 
the semiclassical approximation of the quantum Baker map ^ . Adding to the 
action an integer valued function I(x) does not affect any semiclassical calculation 
(see ||). The choice I(x) — [2x] renders the action invariant under the symmetries 
of the mapping. Inserting an overall minus (a matter of convention) gives: s(x, x') = 
s(x) = (x — l)[2x]. We use the symbol s to denote the actions because the phase- 
space area is \\A\\ = 1, so the actions are properly normalized. Applying the above 
to the i'th periodic point (or segment) of a periodic orbit v of length n, yields 
s(xi,Xi + i) = s(xi) — (xi — l)[2xi\ = (xi — l)vi. The total action of the periodic 
orbit v is the sum of its segment actions 

n n 
i=l i=l 

The action thus defined is invariant under space reflection R, time reversal T and 
the baker transformation B which is only a cyclic permutation of (yx....v n ) i.e. 

s n {v) = s n {R{ V )) = a„(f(i/)) = s n {B{v)) . (24) 

Since the periodic orbits v 1 R{v) 1 T(v), RT(y) are not identical in general, one 
expects a maximal symmetry degeneracy of 4. 

The action can also be expressed in a matrix form Jl9[] which will be useful later. 
Writing (v\ = Vec(v\....v n ) and \u) — Vec{v 1 ....v n ) T then the action is 

s n {v) = {v\6\v)-{v\v) , (25) 

where O is a n x n cyclic matrix, with matrix elements: 

— 1 — 2 (i-j+ n -i) ro odn (26) 
v 2 n 1 
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This expression of the actions clearly shows that their values are restricted to integer 
multiples of 2~ n . This is the minimum separation in the action spectrum. 

The Baker map is uniformly hyperbolic. The stability eigenvalues of all the n- 
periodic orbits are the same, namely A± = 2 ±n , the corresponding eigenvectors are 
parallel to the (x, y) axes, and the Maslov indices are all null. These features bring 
a large degree of simplification which is to be incorporated in the definition of the 
action density for the Baker map. 



3 The action density 

The action density (||) for the Baker map takes the form 

dd(s; n) = 1 6 ( s - «(")) > ( 2? ) 

2 2—2 2 * ' 

where the summation extends over all the different vectors v. Note that in this 
way of writing, repetitions are properly weighted. Since the weights in (|27j) are all 
positive, it is convenient to normalize the density to unit integral. We shall denote 
the normalized density by P n (s), with J P n (s)ds = 1 (the subscript cl is dropped 
since we shall deal exclusively with the classical spectrum). 

P ^ = TT E -^( s " s a) (28) 

where P n stands for the set of distinct n-periodic orbits and r a is the repetition num- 
ber. The normalization factor N n , the number of n-periodic orbits, tends to — in 
the large n limit. An alternative expression is obtained by lumping together all the 
orbits which have the same action and the set V n consists of single representatives 
from each degeneracy set, and g a is the corresponding degeneracy, 

P ^ = W E -*(*-*-) . (29) 
N n r a 

As defined in (||||), P„(s) is a distribution. For some purposes, it is advantageous 
to use a smooth version obtained by convoluting (^8 29) with a narrow window 
function. The resulting smooth function will also be denoted as P n (s), and it is 
shown in figure |. for n = 21. 

The simplest (but wrong) estimate for the action density is obtained by assuming 
that the segment actions s(xj) = (x; — ( p3[ ) are independent random variables 
which arc uniformly distributed on the interval [—^,0]. This leads to a Gaussian 
distribution in the limit of large n, which has the same mean as the action distri- 
bution but otherwise is quite different from it. Figure |]. compares the distribution 
computed numerically for a random choice of segments (n = 19), with the actual 
action distribution (periodic orbits). 

We shall start the discussion of the action density by reviewing some general 
properties. In the next subsection, we shall investigate the substructures in the spec- 
trum, and show that it partitions naturally to families which can be characterized 
in terms of symbolic codes. 



3.1 General properties 

The action distributions P n (s) is characterized by various scales whose dependence 
on n will be summarized below. 
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1.5 -1 -0.5 



Figure 3: Spectrum of actions, n=21 



I I random segments 
I I periodic orbits 



Figure 4: The density of randomly generated actions for n = 19, compared with 
the true density. 



Degeneracy in the action spectrum is mainly due to the symmetries (P4|). As n 
increases, a larger fraction of the orbits are not self symmetric, and therefore the 
mean degeneracy is (g) = 4. Figure (||) shows the degeneracy distributions for 
n=16,17. One can see that for these values of n the degeneracy 2 and 1 still 
appear with appreciable frequency, and higher degeneracies whose origin is number 
theoretical are also possible. 

The lowest scale is 2~™ since the actions are integer multiples of this interval. 
The largest scale is provided by the interval /„ which supports the action distri- 
bution. It follows directly from (|25|) that the maximum value of s n is s™ ax = 
while the minimum is given by the action associated with the periodic orbit v = 
(010101 ...). 



s„ = - 



nii 1 2" + l 

6 3n 2*»-l 



for n even 



for n odd 



Hence, to leading order /„ 



This interval accommodates N n 



2" 



(30) 



periodic 



orbits, which are m 4 times degenerate. It follows that 
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Figure 5: The distribution of degeneracies for n = 16, 17. 



The mean spacing is 



As {n) 



2n2 
32™' 



(31) 



This estimate reproduces very well the numerical simulations. The above result 
indicates that not all the integer multiples of 2~™ on the action line are populated. 
Rather, the gaps are of order n 2 . They cannot be seen in figure (||), because the 
bin size used is too coarse. 

To characterize the large scale features of the action distribution we quote ex- 
plicitly its first 3 m oments in the limit of large n (the proof of these relations is 
given in section ( |3.3| ) below) : 



(s r , 



var(s ra ) 



n 

192 

n 

768 



(32) 



To check the limiting action distribution, it is appropriate to examine its depen- 
dence on the scaled action: 

s — (s n ) 

S = / t \ (33) 
v /var(s„) 

Sano H has recently shown that the scaled action density becomes Gaussian in the 
limit of large n. Our numerical computations confirm that 



lim 



(Sn ~ (S„)) 



n^oo (Z-l)!!n z / 2 



= 1 



lim 

n — *oo 



(Sn - (Sn)) 1 ) 

^75 







for I even 



for I odd 



(34) 



Figure (||) compares the scaled distributions for n = 30, 100, (computed for the r = 5 
level of the ising model of section (3.3). One can clearly see that the distribution for 
the larger n gets more symmetric. The fact that the distribution of scaled actions 
tends to a Gaussian does not imply that pair correlations do not exist. 



3.2 Families of actions and symbolic codes 

Symbolic codes are naturally associated with a partition of phase-space. (See e.g., 
|23J). The actions, being functions of phase-space points are expressed in terms of 
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Figure 6: Scaled action densities for n = 30, 100. 



the codes (|2^) and therefore their classification into families which share certain 
properties is conveniently carried out in terms of their codes. 

As was shown in section |2~2] the binary code {0, 1} is based on the partition of 
phase space by a vertical line at x — \ into the two rectangles which are associated 
to the codes 0, and 1. Only the points of the trajectory which fall in the second 
rectangle (code 1) contribute to the action, and the increment to the total action 
per point is Xi — 1. Denote by p the number of times the periodic orbit visits the 
second interval, or equivalently, the number of z/j — 1 in the code. The set of all the 
trajectories which have the same p will be referred to as a "p-family" . The action of 
trajectories within a p-family cluster in substructures which are illustrated in figure 

in terms of the densities P n , P (s), for n = 23 and a few values of p. The P n ,p(s) 
are narrower than P n (s), and one can show that they are centered about — ~ E^Hl^E) t 
Scaling to unit variance and shifting to their mean s value (figure ^.), the densities 
Pn,p(s) collapse to a single function which is very similar to the scaled total density. 
Thus, the families enable a study of the spectrum of actions with a finer resolution, 
with an approximate scaling similarity of their densities. 



n = 23 p =[ 6, 7, 8, 9, 10, 11] 




Figure 7: distributions of actions according to p-families. 

The partition of the actions into p-families is just the first level in a systematic 
procedure which enables the sorting of the actions according to their codes. The 
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Level 1 n=23 
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8 


-e- p = 


9 


-»- p = 


10 


-B- p = 


1 1 




2 4 

S 



Figure 8: Scaling of the p-family distributions. 

higher the level, the more refined are the resulting families. The higher level families 
still maintain the scaling similarities which was discussed for the p-families. 

The r'th level families, (with r < n) are obtained by a partition of phase space 
into 2 r identical rectangles of unit hight, and with ^f- < x < for all 1 < j < 2 r . 
The family consists of all the rt-periodic orbits which go through each of the 2 r_1 
rectangles in the rightmost half of the unit square the same number of times. Given 

(r) 

a rt-periodic orbit, we denote by p\ , the number of times it goes through the Z'th 
rectangle with \ + < x < \ + and 1 < I < 2 r_1 . A family of the r'th level 
is characterized by the set of 2 r_1 numbers 



P W 



( (r) (r) 
[Pi ,P 2 > 



(r) 



(35) 



(r) 

One can get the P code by direct inspection of the binary code v. p\ ' is the 
number of times the r digits string, beginning by 1 and terminating by j— 1 written 
in a binary basis, occurs in v i.e 



j—l in binary 



/,/ iU} :: ■ [ ()()■■ -111 ■ : mV 



(36) 



r digits 



One can easily check that the p families defined at the beginning of this sections 
correspond to r = 1, with p^' = p. The r and the r + 1 partitions are related by 

(r) (r+l) . (r+1) 
Pj =P2j +P2j-1 

Also, 

Mr : 2^ Pj -Pi = P ■ 
The mean action for a P( r ' family can be written as 



ST (r) ( 1 , 2 j - 1 



2 r+l 



Pi 1 



p E,-=i j p 



(r) 
j 



£.7=1 P. 



(r) 



(37) 



Actions which belong to the same family cluster about this mean value within an 
interval of order . 
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The partition into families break the R symmetry in the sense that symmetry 
conjugate pairs of periodic orbits occur in different families. This is clear from the 
observation that the number of 1 in the two conjugate codes are p and n — p. For 
r > 3, periodic orbits which are related by T can be assigned to different families. 
(Example: the periodic orbits which are represented by the binary codes (001011) 
and (001101) are T conjugate but they appear in the r — 4 families (01011000) and 
(01100100), respectively). To avoid this problem, one should partition the actions 
which correspond to the set V n , where only one representative of each degeneracy 
class is included (^9|). 

An alternative partitioning of the actions in families can be defined in the fol- 
lowing way. Given a binary code v, let q\ be the number of 1 in v, let q 2 be the 
number of sequences {11} in v, and in general, q r is the number of sequences of 
r consecutive 1 in v. (Clearly qi — p). The r'th level family consists of all the 
n-periodic orbits which have the same QW code, 

Q M = (<Zi r U r \---^ r) ) • (38) 

This partitioning has the advantage that T-conjugate orbits are always in the same 
family, (however, i?-conjugate orbits are in different families). We found this parti- 
tioning more convenient for numerical simulations. The two methods of partitioning, 
at the same r level, have in common the same resolution 2~ r . 

The partition of the actions of n-periodic orbits into families provides a very 
useful tool for the analysis of the spectral correlations, especially because the pa- 
rameter r which is at our disposal, determines the resolution at which we wish to 
interrogate the spectrum. Given r, most pairs of actions at a distance less than 
2~ r belong to the same family, and hence, intra family investigations are sufficient 
for the study of smaller action differences. However, for larger scales, all actions 
in a family can be lumped together, and the large range correlations are expressed 
through the study of inter family correlations. We shall take both approaches when 
we discuss the two-point correlations in the action spectra. 



3.3 An Ising model approach 

In the present section we shall introduce a method to investigate the action spectra 
on scales which are larger than a given level of resolution. For this purpose, it is 
convenient to study the Fourier transform of the action distribution, P n (k), 

p nik) = ^E e4fcs " M = ^£ e KE:> o ^-i:;^) (39 ) 

M M 

where the sum is over all possible configurations of the code v = (vx, v%, ■■■.v n ), 
Vi G {0, 1}. The action is given explicitly by 

S K) = 1 _ X 2 _ n +\ v i v i+i + \ v i v i+* + -" ( 40 ) 



Formally ( p9[ ) is the partition sum of a one dimensional Ising model on a circular 
lattice with exponentially decreasing interactions and an imaginary temperature. 
It can be evaluated using the standard transfer matrix method. Using this method 
we can approximate P n (fc) by truncating the interaction at any desired range r, 
(1 < r < n — 1). Thus, r = 1 is the nearest neighbors approximation, r = 2 is 
the next to nearest neighbors approximation etc... . When r = n — 1 we regain 
the full range of interactions. The original effective hamiltonian is invariant under 
space reflection R and time reversal T. The invariance under T is due to the cyclic 
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property of O and hence, is maintained for any truncated version. This is not the 
case however for the space reflection R symmetry. In order to preserve R invariance 
one has to rescale the strength of the term . v 2 which was introduced for R 
invariance of the full action (|2^). For a given range r of the approximation, the 
action takes the form 



1 



1 



n 



£=0 i=l i=l 



(41) 



The symmetry restoring coefficient 7«(r) is calculated by demanding that for 
the r tfl order approximation the action be R symmetric i.e. s n ( r ) 
and is given by 



7: 



(42) 



One recovers 7 ra (0 = 1 for the full range interactions r = n — 1. 

This completes the definition of the r'th level approximants, and because of 
the exponential decreasing strength of the interactions we can expect that any 
quantity computed at the r'th level, will converge exponentially fast to its full 
range interaction value. 

Note: for the rest of this section we omit the factor x J^- n ~ 1 which multiplies 

We express P n (k), for a given range r as: 

Pi r) (k) = ^tr{[T« (*)]» } , (43) 
where the transfer matrix T^ r '(k) is given by 



/ 1 





T {r \k) 



\ 





[3 2 "- 1 [3 2r - 2 
/3 2 "- 3 /3 2 "^ 





1 



\ 





1) 



(44) 



Since T depends on powers of j3, tr j[T( r )]™} is a polynomial in (3 with real coeffi- 
cients Aj 

N(n,r) N(n,r) 



ltr{[T(<-)(fc)]*}= J2 A ^= E A 



e~ lK — (45) 
where N(n,r) is the degree of the polynomial. Transforming back to P n (s) yields 



N(n,r) 



2 r+l- 



(46) 
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and 



N(n,r) 

E 4 } 

3=0 



(47) 



The r'th approximant to the action spectrum consists of N(n, r) equally spaced 
actions Sj — — , each weighted by a normalized weight (or probability) . The 



Aj's contain, together with the spacing Ss(n,r) = ^rrr > all the information about 
the statistical properties of the spectrum at the resolution ^ftt • Increasing the 
range r to r + 1 results in approximately doubling the number of actions, in addition 
to distributing them on a lattice with half the spacing. In the limit r — ■> n — 1 (full 
range interactions) the lattice spacing becomes 4r (or 2 „ 1 _ 1 taking into account 
the factor 



1-2- 



that was neglected). This observation enables us to connect the 
present approach with the partitioning of the actions to families according to their 
codes. The truncation of the interaction at the r'th level is approximately equivalent 
to replacing the actions of all the members of a family by their average value given 
by (|37l). Thus, the coefficients 2 n A^p approximate the sum of the cardinalities of 

all the families whose average action is Sj = — -^+r ■ 

The moments of the distribution are computed using 



'n.(r) 



d m Pi r) (k) 



2 r+l 



dk m 

E 

J=0 



k=0 



(48) 



The two lowest level approximants can be solved analytically. 
r = 1 : The transfer matrix ([l4|) for this "nearest neighbors interaction" approxi- 
mation is 



Its eigenvalues are A± = 1 ± \fj3 which yields 



= e 



(49) 



Pi^Hk) = -L[(i + ^)" + (i-^)' 



• E 

m— 0:everi 



and, 



with 



P, 



2 

2" 



(50) 



(51) 



s e { , 



12 3 1 rn 

"4 ' ~4 ' ~i ' ••• ' ~4 [2 

This distribution limits to a Gaussian distribution for large n. The mean and the 
variance are given by 

/ \ n 1 \ n 

(s„, ( r=i)) = -jg ; var(s„ !(r=1) ) = — 

r — 2 : The transfer matrix in this "next to nearest neighbors interaction" approx- 



imation is 



T 



0=2) 



(fc) 





1 


1 





°\ 










1 


1 




/3 3 


/3 2 








I 








p 


1 1 



; a 



(52) 
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The eigenvalues of T( r_2 ) are 



Ai± = \ (l-/3± VAT) ; A 2± = i (l + /3±VA^ 
with Ai = 1 + 2/3 - 3/3 2 ; A 2 = 1 - 2(3 + 5/3 2 . This gives: 



(53) 



2 2 



n n—m ~ ~ " 



n W n - m \ / t W f _ " 
to 



m—0;even j—0 u—0 t—0 

with 

C jtu = 2*{3"(-iy +u + 5 u (-l)*} . 
The action distribution is 

P(- 2 >( S ) = 



(54) 



2 2 



with 



n n—m ~ p 

m—0;eveji j—0 u—0 



rn 



-8s - 2u - j 



Cj U (s) 



C ju (s) = 2- 8s - 2u - : > {T(-iy +u + 5 u (-l)- Ss - : >} . 
The action takes the values 



12 3 n 

S ^ \ > ~o ' _ 77 > _ 77 ' ••• ' _ 77 



and the summation coefficients u, j, to must fulfill the conditions 



2u 



+ j < -8s < u + j + m/2 



to keep the binomial coefficients well defined. The r — 2 distribution is not Gaussian, 
its mean, variance and third moment can be computed by a straight forward but 
cumbersome calculation. 



71 i Tt I 

(s„, (r=2) ) = - — (1 + -) ; var(a„ i(r=2) ) = — (1 + -) 



(( s n,(r=2) - ( s n,(r=2))) 



3n 
4096 



We were unable to compute analytically the action distributions at higher levels, 
since this involves finding the eigenvalues of the 2 r dimensional transfer matrices. 
However, the computation can be performed using computer codes which perform 
algebraic manipulations. For any r and f3 one can compute tr |[T( r )(fc)]"} by raising 
to the power n and taking the trace. Using Newton's identities, one can 
compute the coefficients of the characteristic polynomial of T 1 - 7 ") from the traces of 
its lower powers, tr {[T^ r ^(A:)] z } , I = 1, • • • ,2 r . The traces of higher powers can 
then be expressed recursively in terms of these coefficients. This way we were able 
to perform computations up to the level r — 9, without reaching the limits of our 
computer resources. 

At the beginning of the chapter we argued that due to the exponential decrease of 
the interactions, computing any quantity at a given level r, converges exponentially 
fast to the full r = n — 1 calculation. Thus, it can be expected that any quantity 
of interest will be given by a geometric series in 2~ r . Knowing the first two levels 
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(r = 1, 2) may allow us to extrapolate, and obtain the leading terms for r = n — 1. 
Applying this strategy for the first three moments, we obtained the expression 




These expressions were checked by comparing them to the computer aided, large r 
calculations. 

We shall return to the results of the Ising model when we discuss the pair 
correlations in the action spectrum which is the subject of the next section. 



4 Pair correlations in the actions spectrum 

So far we discussed the distribution of actions. Now we turn to study their pair 
correlations which is the main issue of the present work. We examine whether 
the classical actions spectrum exhibits correlations, and whether these correlations 
conform with the semiclassical theory which was presented in section (|l|) . 

There are two distinct length scales in any discrete and finite spectrum: the 
width or support of the distribution on the largest scale, and the nearest neighbors 
spacing on the smallest scale. The fact that the support is finite induces a trivial 
correlation: the probability to find two points separated by a distance larger than 
the range of the support is zero. The nearest neighbors scale is of great interest, 
in particular for "rigid" spectra such as the quantum spectra of classically chaotic 
systems. They exhibit level repulsion as predicted by RMT i.e. the probability to 
find two nearest neighbors at distance 6 vanishes as S — > 0. We have already hinted 
in section (|l]) that the pair correlations in the action spectrum is expected on a 
scale which is much larger that the mean spacing. 

We shall perform the analysis of pair correlations in several steps. In the first 
we shall examine the correlations on the scale of the mean level spacing. We shall 
show that the m'th nearest neighbor spacing distribution are consistent with the 
corresponding Poisson distributions, as long as m is smaller than the correlation 
length measured in units of the mean spacing. We shall then study the formfactor 
of the pair correlation function, and compute it in various ways. The main result 
of this investigation is that the correlations exist, and their scaling with n agrees 
with the semiclassical expectations. 

4.1 The m'th- neighbors spacings distributions 

The m'th- neighbors spacing distribution p(")(cr;m), is a straight-forward general- 
ization of the widely used nearest neighbors spacing distribution. To define this 
distribution, the action spectrum is ordered so that Sj < Sj+i, degeneracy sets 
are represented by a single value and the spectrum is normalized by the nearest 
neighbor spacings 

Then, p^ (c; m) is the probability that j^^y takes the value a. 
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Figure 9: The first, second and third neighbors distribution (n = 19) 



As a reference spectrum we generated numerically a Poissonian spectrum. It is 
randomly chosen from a Gaussian probability density which has the same mean, 
variance and degeneracy structure as the action spectrum for the n value of interest. 
In the upper frame of figure ||. the nearest neighbors distribution, p i - n \a- 1 m = 1), 
for n — 19 is plotted together with the random (Poissonian) nearest neighbors 
distribution. The similarity between the two distribution, which persists over 6 
orders of magnitude, provides strong evidence in favor of the claim that on the 
mean spacing scale, the action spectrum is statistically random. 

This finding is further corroborated in the lower frames of figure |9[ were the 
spacing distributions for m = 2, 3 are compared with the corresponding random 
distributions. 

To go even further in m, we studied numerically the summed distributions 

M 

p { m\°)= $> ( "V;™) , (59) 

m— 1 

which approximates the two point correlation function R2- n ((r) — Poo(c) in the 
range < a < M. 



t + periodic orbits, n=1 9 

& random actions 
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Figure 10: Integrated distributions with (M = 20, 50, 100) , (n = 19) 
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Figure [To], shows p$ '(a) for M = 20,50,100, (n = 19) together with the 
corresponding random distributions. These values M should be compared with the 
correlation length expected to be of order 

Ucorr " (A^O) = 3 ~ («9Qforn = 19). (60) 

This numerical investigation clearly demonstrates that the systematic deviation 
between the random and the actual distributions increases as M approaches the 
correlation length. 



4.2 The classical spectral formfactor 

The purpose of the present section is to study the formfactor of the classical action 
spectrum (|l^). We start by investigating alternative methods for averaging the 
formfactor, the need for which, and the difficulties involved, were discussed at length 
in section (Q). We shall study the function /C n (fc) (|l7|), and in particular test whether 
the action spectrum of the Baker map, satisfies either of the equivalent relations 



IC n (k) 



or , rfCr, 



2irn 

T 



Krmt (t ) . 



(61) 



Since the Baker map is invariant under time reversal, the appropriate RMT expres- 
sion is the one for the Circular Orthogonal Ensemble (COE) Kcoe(2t), where the 
factor 2 in the argument is due to the invariance of the Baker map under the R 
symmetry. The explicit expression for Kqoe (t) is given by 



( 2r-Tln(l + 2r) for r < 1 
' 2-Tlnff^i) forr>l 



(62) 



Before applying any averaging, the function r/C„ 
ations which are shown in figure 
result. The large fluctuations make the comparison quite meaningless, and the 
figure is shown to emphasize the need of averaging. 



displays strong fluctu- 
for n = 15 together with the expected COE 




Figure 11: The function r/C„ (^p) (n=15) ([n]) and the RMT prediction 
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Using the running average: 



f(y)&y , 



(63) 



the large fluctuations are reduced, and the curve labelled "periodic orbits" in figure 
fl2| . is the running average obtained by computing the formfactor from the set of 
n = 17-periodic orbits. The agreement with the corresponding COE curve persists 
up to r w 0.5. The line marked "diag" is the curve obtained by assuming that the 
actions are not correlated, and it agrees very well with the line marked as "rand" 
which was computed for the random set of actions (see section (}4.l|)). The differ- 
ence between the data and the random curves is a clear indication of the presence 
of correlations, whose similarity to the predicted COE result goes beyond the lead- 
ing "diagonal" approximation. The running average procedure is not satisfactory, 
because it is practical only for low values of n. As n increases, the number of pe- 
riodic orbits proliferates exponentially, and the fluctuations in the formfactor grow 
as rapidly. Moreover, this method does not help to unravel the dynamical origin of 
the correlations. 



1.5r 




Figure 12: A running average @ of r/C„ (^) (n=17) 



A systematic averaging procedure, which can be applied for any n, uses the 
concept of families which was introduced in section (|]^) . The spectrum of actions of 
n-periodic orbits is partitioned to families at the r'th level with labels P( r ) (see (Sjf)). 
As long as r is sufficiently small such that 2 _r is larger than the correlation length, 
only intra family correlations are important. The formfactor can be approximated 
by an incoherent sum over the formfactors which pertain separately to orbits within 
one family. 



1 



E 

PW 



E • 



(64) 



Each family of actions is supported on an interval of size which decreases with r as 
2~ r . Hence from the requirement kAs > 2tt we deduce N > 2 r , and the range of 
values of t which can be described by this method is bounded from above by w . 
Thus by increasing r, we gain a higher level of smoothing, but we lose on the range 
of r where this method can be used. 

The resulting formfactors for n = 15 and r = 1,2,3 are shown in figure |l3|. It 
shows that the smoothing gets more effective as r increases, without appreciable 
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Figure 13: family averaged r/C„ (^) for r = 1, 2, 3 , n = 15 



loss of correlations in the interval < r < 1. This is consistent with the expectation 
that the spectral correlation length is which is smaller than the family separation 
of w 2~ r for the values r = 1,2,3. A closer comparison of the r = 2 and r = 3 
curves near r = 1 shows that the deviations from the RMT prediction starts earlier 
for the r = 3 data. To investigate this trend further, we use the Q( r ) partitioning 



(38) which is more convenient from the numerical point of view. It allows us to 
extend the level further, and the data for r = 4, 5 is shown in figure [l4|. Since now 
the family range is smaller than ^ , the correlations can be studied only on a lower 
range of r value. This, and similar other numerical investigations provide strong 
evidence in support of the i dependence of the correlation length, in agreement 
with the semiclassical expectations. 




r=5 



COE 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

T 

Figure 14: Q( r ) family averaged r/C n for r = 4, 5 and n = 15 

As was commented previously, the partition into families breaks pairs of R 
conjugate orbits. Hence, the RMT expression which is used in figures , |l4|. 
is i^cOE(T) and not ifcOE(2r) which is used for the comparison with the full data 
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set. 

The Ising model approach (3.3) offers a complementary view of the formfactor. 
As was explained in (3.3), the truncation at the r'th level does not allow us to 
distinguish between actions which are within 2 _r , and counts them as if they are 
degenerate. Thus, the only way by which we can get spectral information on the 
scale of i in this model is to use r such that \» 2~ r . The spectral correlations 
are now studied as correlations between degeneracy groups, rather than within a 
group. Figure [l5]. shows the average of a few form factors, corresponding to values 
of n from n = 17 to n — 22, given by the model with r = 9 which satisfies the 
criterion above. The average over n reduces the fluctuations of the form factor so 
it is possible to compare it directly to the RAIT. 




Figure 15: rJC n , from the Ising model with r = 9, averaged over 17 < n < 22. 



The fact that the action spectrum is resolved with an accuracy of 2~ r sets a 
lower limit to the values of t which are accessible by the model, which is r m i„ w Tpr . 
This is why the low r values in figure Flq. are absent. 



5 Discussion 

The semiclassical predictions which were presented in section (|l|) were checked in 
detail and were found to be well satisfied for the Baker map. They consist of the 
following main points: 

• The action spectrum show pair correlations which are universal and consistent 
with the predictions of the semiclassical theory and RMT. 

• The correlation length is of order i which exceeds the mean spacing by a factor 
which grows exponentially with n. The spectrum of actions is Poissonian on 
smaller scales. 

• The correlations can be associated with families of periodic orbits which have 
a similar dynamical structure which can be associated systematically with 
their symbolic codes. 

Various aspects of these points were confirmed for other systems [g, [| 0, f| . 
The results of f| § are unique, because the phase space of the mapping considered 
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is four dimensional, and the semiclassical theory predicts that the action correlation 
range scales with nas 4- This scaling was confirmed. 

Richter and Sieber [jllj pointed out the pairs of periodic orbits which provide the 
correlations which are necessary for the second order term in the small r expansion 
of the formfactor. These orbits spend about half of their time quite close to each 
other, and the rest of the time, they go along the time reversed part of their partner's 
orbit. This way both orbits visit the same parts of phase space, but the order 
by which they do it is different. Similar pairs were also proved to provide the 
answer in the case of quantum graphs fl3|| . The families which we identified here 
as storing the correlations between the actions are of the same nature, but of a 
more general character. They do not rely exclusively on pairs conjugated partially 
by time reversal (as explained above), but by the requirement that the correlated 
members spend the same amount of time in the same sub domains of phase space. 
This definition of families is also consistent with the work of [Q , where the periodic 
orbits in a family follow the same segments of phase space which were related to 
each other by symmetries other than time reversal. Because of the strong link 
between phase space partition and symbolic codes, it was possible to characterize 
the families by a common code. Unfortunately, even with this tool, we were not 
able to derive the classical correlation function using classical arguments only. This 
remains an enigma which should be addressed. At the same time, one should also 
be able to show that the spectrum of actions on smaller scales is Poissonian. This 
is assumed, but not discussed in [jllj and in p3[ . 

In the present article we focussed our attention on discrete maps and their peri- 
odic orbits. However, most system of interest are Hamiltonian flows, where the time 
is a continuous variable, and the quantum spectrum is on the real line and not on 
the unit circle. This does not pose any essential problem, since several quantization 
techniques make use of an auxiliary map to derive the quantum energies, and it 
was shown that for chaotic systems, the spectral statistics of the energies and of the 
eigenphases of the auxiliary map are the same in the semiclassical limit |^5|, |2(| || . 
A more direct approach was recently introduced in [^7| , where the quantization of 
the hamiltonian flow is carried out in terms of a quantum map which evolves the 
system along a sequence of equally spaced times t n — n At. The semiclassical 
expression for the spectral formfactor is analogous to (0), and the periodic orbits 
have periods which are integer multiples of At and their energies are restricted to 
a well defined energy interval. 

Integrable dynamics lead to Poissonian spectra [^8| . In the present context this 
implies that the actions are Poissonian too JtJ . Recently it was shown in |^9| that in 
order to account for finer spectral correlations which are due to the spectrum being 
pure point, finer correlations must exist, but this discussion exceeds the scope of 
the present paper. 
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